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Abstract. We study experimentally the thermal fluctuations of energy input and 
dissipation in a harmonic oscillator driven out of equilibrium, and search for 
Fluctuation Relations. We study transient evolution from the equilibrium state, 
together with non equilibrium steady states. Fluctuations Relations are obtained 
experimentally for both the work and the heat, for the stationary and transient 
evolutions. A Stationary State Fluctuation Theorem is verified for the two time 
prescriptions of the torque. But a Transient Fluctuation Theorem is satisfied for the 
work given to the system but not for the heat dissipated by the system in the case of 
linear forcing. Experimental observations on the statistical and dynamical properties of 
the fluctuation of the angle, we derive analytical expressions for the probability density 
function of the work and the heat. We obtain for the first time an analytic expression 
of the probability density function of the heat. Agreement between experiments and 
our modeling is excellent. 
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1. Introduction 

Nanotechnology, as well as biology, biophysics and chemistry are using or studying 
setups and objects which are smaller and smaller. In these systems, one is usually 
interested in mean values, but thermal fluctuations play an important role because 
their amplitude are often comparable to the mean values. This is for example the case 
for quantities such as the energy injected in the system or the energy dissipated by the 
system. These fluctuations can lead to unexpected and undesired effects : for instance, 
the instantenous energy transfer can be reversed by a large fluctuation, leading energy 
to flow from a cold source to a hot one. These events, although rare, are quantitatively 
studied by the recent Fluctuations Theorems (FTs). These theorems give Fluctuation 
Relations (FRs) that quantify the probability of these rare events in systems which can 
be arbitrarily far from equilibrium. FTs have been demonstrated in both deterministic 
systems [TJ [2] and stochastic dynamics [3], HI El [7J [8]. Experiments searching for 
FRs have been performed in dynamical systems P, [101 E] , but interpretations are very 
difficult because a quantitative comparison with theoretical prediction is impossible. 
Other experiments have been performed in stochastic systems described by a first order 
Langevin equation: a Brownian particle in a moving optical trap p2| and an out-of- 
equilibrium electrical circuit [13] in which existing theoretical predictions [6l [7J were 
verified. Interesting comments on the Langevin equation can be found in [T3] , 

In the present article, we study a thermostated harmonic oscillator described by a 
second order Langevin equation. We experimentally search FRs for the work done by 
an external operator and for the heat dissipated by the system, and present analytical 
derivations of FTs based on experimental observations. 

This paper is organized as follows. In sectionEJ we present the experimental system, 
write its energy balance to define the work given to the system together with the heat 
dissipated. We then introduce the Fluctuation Relations (FRs) and the Fluctuation 
Theorems (FTs). In sections |3l H] and [5l we present experimental results on the 
fluctuations of first the work and then the heat. A short discussion on experimental 
results in given in EJ Then, in sections [7J and El we present some analytical derivations 
of FTs based on hypothesis inspired by experimental observations. We compare these 
analytical predictions to the experimental observations and finally conclude in section 
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Viscous fluid 




Figure 1. a) The torsion pendulum, b) The magnetostatic forcing. 



2. System description 

2.1. The harmonic oscillator 

Our system is a harmonic oscillator and we measure the non-equilibrium fluctuations 
of its position degree of freedom. The oscillator is damped due to the viscosity of 
a surrounding fluid that acts as a thermal bath at temperature T. Our oscillator, 
depicted in Fig. is a torsion pendulum composed of a brass wire (length 10 mm, 
width 0.5 mm, thickness 50 /im) and a glass mirror glued in the middle of this 
wire (length 2 mm, width 8 mm, thickness 1 mm). The elastic torsional stiffness of 
the wire is C = 4.65 • 10 -4 N.m.rad -1 . It is enclosed in a cell filled by a water- 
glycerol mixture at 60% concentration. The system is a harmonic oscillator with 
resonant frequency f Q = \JC / I c r/ [2-k) = uq/{2tx) = 217 Hz and a relaxation time 
r a = 2I c fi/v = 1/a = 9.5 ms. / c fj is the total moment of inertia of the displaced 
masses (i.e. the mirror and the mass of displaced fluid) [15]. The damping has two 
contributions : the viscous damping v of the surrounding fluid and the viscoelasticity 
of the brass wire which can be neglected here. 

The angular displacement of the pendulum 9 is measured by a differential 
interferometer [16], [T7J [18]. The measurement noise is two orders of magnitude smaller 
than thermal fluctuations of the pendulum. 9{t) is acquired with a resolution of 24 
bits at a sampling rate of 8192 Hz, which is about 40 times f . We drive the system 
out-of-equilibrium by forcing it with an external torque M by means of a small electric 
current J flowing in a coil glued behind the mirror (Fig. [TJo). The coil is inside a static 
magnetic field. The displacements of the coil and therefore the angular displacements 
of the mirror are much smaller than the spatial scale of inhomogeneity of the magnetic 
field. So the torque is proportional to the injected current : M = A.J ; the slope A 
depends on the geometry of the system. 
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The angular displacement 9 of this harmonic oscillator is very well described by a 
second order Langevin equation: 

d 2 9 d6 



I eB -^ + u — + C6^M + y/2k B Twn, (1) 

where r\ is the thermal noise, delta-correlated in time of variance 1 and k B the Boltzmann 
constant and T the temperature of the system which is the one of the surrounding fluid. 
The Fluctuation Dissipation Theorem (FDT) gives a relation between the amplitude 
of the thermal angular fluctuations of the oscillator at equilibrium and its response 
function. For a harmonic oscillator, the equilibrium thermal fluctuation power spectral 
density (psd) is: 

/mi2\ ^bT a 4k B Tu 

m > = _ [mx = ^————^ (2) 

where x = f = Aj. Using FDT (Eq. 0, we measure the coefficient A and test the 
calibration accuracy of the apparatus which is better than 3%. More details on the 
set-up can be found in [161 Ej- 

2.2. Energy balance 

When the system is driven out of equilibrium using a deterministic torque, it receives 
some work and a fraction of this energy is dissipated into the heat bath. Multiplying 
Eq. (fTJ) by ^ and integrating between tj and tj + r, we obtain a formulation of the first 
law of thermodynamics between the two states at time tj and tj + r (Eq. ([3])). The 
change in internal energy AU T of the oscillator over a time r, starting at a time ti, is 
written as: 

AU T = U(U + r) - U(U) = Q T + W T (3) 
where W T is the work done on the system over a time r : 

and Q T is the heat given to the system. Equivalently, (— Q T ) is the heat dissipated by 
the system. AU T , W T and Q T are defined as energy in k B T units. The internal energy 
is the sum of the potential energy and the kinetic energy : 

m = kir ^ /cff + ^ ce ^ 2 y ( 5 ) 




The heat transfer Q T is deduced from equation ([3]) ; it has two contributions : 

Q T = AU T - W T 

1 f u+T \d6 /n ]\, 1 f u+T . n d9, ni , 

= -kPrl "U (i ) J it + klr.k " (i V t )dt ' (6) 

The first term corresponds to the opposite of viscous dissipation and is always negative, 
whereas the second term can be interpreted as the work of the thermal noise which have 
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a fluctuating sign. The second law of thermodynamics imposes (—Q T ) to be positive. 
We rescale the work W T (the heat Q T ) by the average work (W T ) (the average heat (Q T )) 
and define: w T = (q T = tqt)- The brackets are ensemble averages. In the present 
article, x T , respectively X T , stands for either w T or q T , respectively W T or Q T . 



2.3. Fluctuation Theorems and Fluctuation relations 

There are two classes of FTs. The Stationary State Fluctuation Theorem (SSFT) 
considers a non-equilibrium steady state. The Transient Fluctuation Theorem (TFT) 
describes transient non-equilibrium states where r measures the time since the system 
left the equilibrium state. A Fluctuation Relation (FR) examines the symmetry of the 
probability density function (PDF) p(x T ) of a quantity x T around ; x T is an average 
value over a time r. It compares the probability to have a positive event (x T = +x) 
versus the probability to have a negative event (x T = —x). We quantify the FR using a 
function 5" (symmetry function) : 

S(x T ) = -^ln( P ^ = +X \ ). (7) 



(X T ) \p(x T = -x) 

The Transient Fluctuation Theorem (TFT) states that the symmetry function is 
linear with x T for any values of the time integration r and the proportionality coefficient 
is equal to 1 for any value of r. 

S(x T ) = x T , Vx T , W. (8) 

Contrary to TFT, the Stationary State Fluctuation Theorem (SSFT) holds only in the 
limit of infinite time (r). 

lim S(x T ) = x T , Wx T . (9) 

The questions we ask are whether fluctuation relations for finite time satisfy the two 
theorems and what are the finite time corrections. In a first time, we test the correction 
to the proportionality between the symmetry function S(x T ) and x T . In the region where 
the symmetry function is linear with x T , we define the slope E x (r) : S(x T ) = T: x (t)x t . 
In a second time we measure finite time corrections to the value S x (r) = 1 which is the 
asymptotic value expected by the two theorems. 



3. Transient non-equilibrium state 

For the Transient Fluctuation Theorem, we choose the torque M(t) depicted in Fig. [2^). 
It is a linear function of time : M(t) = M t/r r with M Q = 11.28 pN.m and 
r r = 0.1 s = 10.52 Tq,. The value of M is chosen such that the mean response of 
the oscillator is of order of the thermal noise, as can be seen in Fig. [2b) where 9(t) is 
plotted during the same time interval of Fig. [2^). The system is at equilibrium at ti — 
(M(ti = 0) = pN.m and M(t) = pN.m Vt < ti). In this section the starting time U 
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Figure 2. a) Typical driving torque applied to the oscillator; b) Response of the 
oscillator to the external torque (gray line) . The dark line represents the mean response 
9{t) to the applied torque M(t). 

of integration of all quantities defined before (W T , AU T and Q T ) is tj = 0. So the work 
is : 

w ^^r[ M KtW- (10 » 

3.1. Average values informations 

In Fig. [3k), we represent the time average ((t~ 1 W t )) of the power injected into the 
system, the internal energy difference (t~ 1 AU t ) and the time average ((r _1 WV)) of the 
power dissipated by the system. (t^Wt) and (t~ 1 AU t ) are linear in r after some short 
relaxation time r a defined in the Langevin equation : for r/r a smaller than 1, some 
oscillations around the linear behavior can be seen. The average value of work (W T ) 
is therefore quadratic in r and is equal to 33 ksT for r = r r . The difference between 
(W T ) and (AU T ) corresponds to the mean value of dissipated heat {— Q T ) (Eq. (jSJ)). As 
can be seen in figJBK), (W T ) is larger than (AU T ) for all times r. The average of the 
dissipated power ((— t~ 1 Q t )) is therefore positive for all times r as expected from the 
second principle. For r larger than several r a , the dissipated power is constant and 
equal to a few ksT per second because (t~ 1 W t ) and (r~ 1 AL r T ) have the same slope 
after some r a . So we have the following behavior : the work done by the external work 
is used by the system to increase its internal energy but a small amount of energy is 
lost at a constant rate by viscous dissipation and exchange with thermostat. 

3.2. Work fluctuations 

The probability density functions (PDFs) p(W T ) of w T is plotted in Fig. [3b) for different 
values of r/r a . Four typical value of r are presented : the first is smaller than the 
relaxation time and the last equals five relaxation times ; the results are the same 
for any value of r. The PDFs of w T are Gaussian for any r. We observe that w T takes 
negative values as long as r is not too large. The probability of having negative values of 
w T decreases when r is increased. From the PDFs, we compute the symmetry functions. 
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Figure 3. TFT. a) Average value of t^Wt (o), t~ 1 MJ t (□) and t _1 Q t (o) plotted 
as a function of r. b) PDFs of w T for various t/t q : 0.31 (o), 1.015 (□), 2.09 (o) and 
4.97 (x). Continuous lines are theoretical predictions with no adjustable parameters. 

c) Corresponding functions S(w T ). The straight continuous line is a line with slope 1. 

d) PDFs of t^ 1 AU t for two values of r/r Q : 4.97 (o) and 8,96 (□). e) Corresponding 
PDFs of q T . Continuous lines are Gaussian fits, f) Corresponding functions S(q T ). 
The straight continuous line is a line with slope 1. 
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They are plotted in Fig. as a function of w T . They all collapse on the same linear 
function of w T for any r, which implies they all have the same slope S w (r). The straight 
line in fig Eb) has slope 1. Within experimental error bars, T< w (t) is equal to 1 for all 
time r. Therefore work fluctuations for a harmonic oscillator under a linear forcing 
satisfy the TFT. We checked that this property is true for other values of M D and r r . 

3. 3. Heat fluctuations 

The PDFs of r~ 1 A[/ T are plotted in Fig. [3U) for two values of r/r a : they are not 
symmetric and have exponential tails. The PDFs of q T can be seen in Fig. [3b) for the 
same values of r/r Q . They are qualitatively different from the ones of the work. We 
have plotted in the same figure the Gaussian fit of the two PDFs of the dissipated heat. 
It is clear that the PDFs of q T are not Gaussian. Extreme events of q T are distributed 
on exponential tails. These tails can be interpreted noticing that Q T = AU T — W r and 
AU T have exponential tails. The variance of the PDFs of q T is also much larger than 
the variance of the PDFs of w T . 

We plot on Fig. [3j) symmetry functions S(q T ) for the same times r/r a . Only the 
behavior of large events can be analyzed here because the variance is much larger than 
the mean a Wr 1. As it can be seen in Fig.EF, S(q T ) is not proportional to q T , therefore 
TFR is not satisfied for finite time. Within experimental resolution, S(q T ) is constant 
for extreme events and equal to 2. This behavior can be interpreted by writing for large 
q T , p(q T ) = A± exp(— a ± |g T |) where a + and a_ are the decrease rate on the exponential 
tails. Each coefficient depends on r. There is a simple expression of S(q T ) for large 
fluctuations: 



In Fig. [3b), the PDFs of q T are symmetric around the mean value for the two values of 
t. It is not the case for small t/t q . Thus we can conclude that a + = a_ and that the 
symmetry function is so equal to the constant : ((Q r )) _1 (ln(v4 + ) — ln(A_)). 

As it can be seen in Fig. [3b), the PDFs become more and more Gaussian when r 
tends to infinity. It is expected that for infinite time, the PDF of q T is a Gaussian. Thus, 
TFT appears to be satisfied experimentally in the limit of infinite r. Our interesting 
finding is that, for Q T TFT if not valid for any times. 

4. Steady state : linear forcing 

4-1. Definition of the work given to the system 

We call a steady state a state in which both forcing and response to the forcing do not 
depend on the initial time ti, but only on r. This implies that (M(tj + r)) i ndependant 
of ti ; and so is (6(U + r)). If the torque drifts along time, the mean of M{U + r) is 
linear with ti + r. Thus we have to change the definition of the work done on the system 




(11) 
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to be in a steady state. This is equivalent to a Galilean change of reference frame. The 
work is now defined as : 



With this definition, the forcing is M(t) — M(tj) and the response to the forcing 
8(t) — 6(ti). When we impose a forcing linear in time (M(t) = M t/r r ), the first 
condition ((M(tj+r)) independent of U) is satisfied. The second ((9(U + t)) independent 
of ti) is also satisfied if tj > 3r a , i.e. after a transient state. Thus the system is on a 
steady state. We remark that, in the transient state, this definition of the work reduces 
to the usual one, because M(ij) = M(t = 0) = pN.m. 

4-2. Work fluctuations 

The average of W T is quadratic in r for any value of r/r a . There are no oscillations 
in time for small r/r a . The PDFs of w T are Gaussian for any value of r/r a (Fig. Wfa). 
Probability of negative values is high and decreases with r, like in the transient case. 
The symmetry functions S(w T ) are again proportional to w T (Fig. [4b ) but the slope T, w 
is not equal to 1 for smaller r and tends to 1 for r ^> r a only, as can be seen in Fig. [4b. 
Thus we obtain a fluctuation relation for the work done on the system in this steady 
state and this relation satisfy the SSFT. The slope at finite time is slightly oscillating 
at a frequency, close to fa. 

4-3. Heat fluctuations 

The heat dissipated during this linear forcing has a behavior very similar to the one 
observed in the transient case (section 13. 3j) . We can so transpose here all what we said 
in section [3731 

5. Steady state : sinusoidal forcing 

We now consider a periodic forcing M(t) = M sm(uidt). This is a very common kind 
of forcing but it has never been studied in this context. Using Fourier transform, any 
periodical forcing can be decomposed in a sum of sinusoidal forcing. We explain here 
the behavior of a single mode. We choose M = 0, 78 pN.m and uj^j (2n) = 64 Hz. This 
torque is plotted in Fig. [5^- The mean of the response to this torque is sinusoidal, with 
the same frequency, as can be seen in Fig. [5b. We studied other frequencies Ud- The 
system is clearly in a steady state. We choose the integration time r to be a multiple 
of the period of the driving (r = 2nn /ojd with n integer). The starting phase t^d is 
averaged over all possible t« in order to increase statistics ; in the remaining of this 
section, we drop the brackets (.) ti - 




(12) 
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Figure 4. SSFT with a ramp forcing, a) PDF of W T for various r/r Q : 0.019 (o), 0.31 
(□), 2.09 (o) and 4.97 (x). b) Corresponding functions S(W T ). c) The slope S w (r) 
of S(W T ) is plotted versus r (□: experimental values; continuous line: theoretical 
prediction eq. (|20p with no adjustable parameters). 




Figure 5. a) Sinusoidal driving torque applied to the oscillator, b) Response of 
the oscillator to this periodic forcing (gray line) ; the dark line represents the mean 
response (6(t)). 
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Figure 6. Sinusoidal forcing, a) PDFs of the work w n integrated over n periods of 
forcing, with n = 7 (o), n = 15 (□), n = 25 (o) and n = 50 (x ). b) The function S(w n ) 
measured at L0d/2ir = 64 Hz is plotted as a function of u>„ for several n: (o)n = 7; 
(□)n = 15 (o)n = 25; (x)n = 50. For these two plots, continuous lines are theoretical 
predictions with no adjustable paramaters (Eq. and Eq. (|.14p ). c) The slopes 

£i/,(n), plotted as a function of n for two different driving frequencies uj c i = 64 Hz 
(□) and 256 Hz (o); continuous lines are theoretical predictions from Eq. (j2Tj) with no 
adjustable parameters. 



5. 1 . Work fluctuations 

The work is written as a function of n, the number of periods of the forcing : 



The PDFs of w n are plotted in Fig. Eh- Work fluctuations are Gaussian for all values of 
n as in previous cases. Thus symmetry functions are again linear in w n (Fig. [6b). The 
slope E w (n) is not equal to 1 for all n but there is a correction at finite time (Fig. [6b). 
Nevertheless, T, w (n) tends to 1 for large n, so SSFT is satisfied. The convergence is very 
slow and we have to wait a large number of periods of forcing for the slope to be 1 (after 
30 periods, the slope is still 0.9). 

This behavior is independent of the amplitude of the forcing M and consequently 
of the mean value of the work (W n ). The system satisfies the SSFT for all forcing 
frequencies Ud but finite time corrections depends on uj d , as can be seen in Fig. [6t. 




(13) 
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n 

Figure 7. Sinusoidal forcing, a) Average value of W n (o) and Q n (□). In the next 
plots, the integration time r is a multiple of the period of forcing, r = 2mr /cod, with 
n = 7 (o), n = 15 (□), n = 25 (o) and n = 50 (x). Continuous lines are theoretical 
predictions with no adjustable parameters, b) PDFs of AU T . c) PDFs of q T . d) 
Symmetry functions S(q T ). e) The slope Ti q (n) of S(q T ) for q T < 1, plotted as a 
function of n (o). The slope S TO (n) of S(w T ) plotted as a function of n(D). Continuous 
line is theoretical prediction. 

5. 2. Heat fluctuations 

We first do some comments on the average values. The average of AU T is obviously 
vanishing because the time r is a multiple of the period of the forcing. (w n ) and 
(q n ) have consequently the same behavior and they are linear in r, as can be seen in 
Fig. [7h) but the PDFs of heat fluctuations q n have exponential tails (Fig. [7b). This can 
be understood noticing that, from Eq. (JSJ), — Q T = W T — AU T and that AU T has an 
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exponential PDF independent of n (Fig. [7b). Therefore, in a first approximation, the 
PDF of q T is a convolution between an exponential distribution (PDF of AU T ) and a 
Gaussian distribution (PDF of w T ). 

Symmetry functions S(q n ) are plotted in Fig. [7H) for different values of n ; three 
different regions appear: 

(I) For large fluctuations q n , S(q n ) equals 2. When r tends to infinity, this region 
spans from q n = 3 to infinity. 

(II) For small fluctuations q n , S(q n ) is a linear function of q n . We then define S g (n) 
as the slope of the function S(q n ), i.e. S(q n ) = S g (n) q n . This slope is plotted in Fig. [Tfe) 
where we see that it tends to 1 when r is increased. So, SSFT holds in this region II 
which spans from q n = up to q n = 1 for large r. 

(III) A smooth connection between the two behaviors. 

We observe that T> w (n) matches experimentally S g (n), for all values of n (Fig. [Tfe). 
So the finite time corrections to the FT for the heat are the same than the ones of FT 
for work : S m (n) = £ 9 (n). 

These regions define the Fluctuation Relation from the heat dissipated by the 
oscillator. The limit for large r of the symmetry function S(q T ) is rather delicate and 
we will discuss it in section 18.21 

6. Discussion and conclusion on experimental results 

In the previous sections, we have presented experimental results on a harmonic oscillator 
driven out of equilibrium by an external deterministic forcing M. We operated with 
two different time-prescriptions : one in which M is a linear function of time, and one 
in which M is a sinusoidal function of time. 

The energy injected into the system is the work W of the torque M. The PDFs of 
the work W are Gaussian whatever the time prescription of M is, and work fluctuations 
satisfy a TFT (M linear in time) and a SSFT (M linear or sinusoidal in time). 

The energy dissipated by the system is represented by the heat Q, and we measured 
it using the first principle of thermodynamics (eq. [6]). Heat probability distributions are 
not Gaussian and are very different from the ones of the work. They nevertheless satisfy 
a SSFT in both the case of a sinusoidal forcing and a linear forcing. But they do not 
satisfy a TFT in the case of a linear forcing, because the symmetry functions are not 
linear for all values of dissipated heat q T . 

In the next two sections, we use some experimental evidences to derive analytical 
expressions of the PDFs of work and the heat exchanged on an arbitrary time interval 
r. We then derive FTs together with their finite time corrections. 

7. Work fluctuations : theoretical predictions 

In this section, we derive the analytical expression of the PDF of the work given to the 
system, and defined as the work of the torque applied to the pendulum, which is either 



14 



linear or sinusoidal in time. Experimentally, we observed that the PDFs are always 
Gaussian, so we restrict our task to deriving expressions for the first two moments of 
the work distribution. 

To do so, we use experimental observations on the fluctuations of the angle 9, 
as described in section 17.11 below. We then compute in section 17.21 the mean and the 
variance of the work W T in the different experimental situations, and then write formally 
the corresponding Fluctuations Relations, from which we obtain analytical expressions 
of the finite time corrections to the Fluctuation Theorems. 

7.1. Angular fluctuations in the presence of forcing. 

We discuss here the angular fluctuations. We decompose the angle 9 into a mean value 
(6) and a fluctuating part 59, writing 9 = (9) + 59. The mean value corresponds to an 
ensemble average. It is obtained experimentally by averaging over realisations of the 
forcing, and it is presented in Fig. [2] and [51 

A first experimental observation is as follows. The measured mean response (9) is 
exactly equal to the solution of the deterministic second order equation obtained when 
removing the noise term (rj = 0) in the Langevin equation ([I]). We checked this from 
our data, and found this way a value of the calibration A (see section 12.11) in perfect 
agreement with the one obtained from the application of the Fluctuation Dissipation 
Theorem. 

A second experimental observation concerns the probability distribution of 59 in 
out-of-equilibrium conditions. We know and observed that at equilibrium, 59 has a 
Gaussian distribution with variance <Jq = ksT/C, and the associated momentum 9 has 
fluctuations 59 which also have a Gaussian distribution, with a variance ksT '/ I e g. We 
observe that the statistical properties of angular fluctuations 59 when a torque M(t) 
linear in time is applied are the same as the statistical properties at equilibrium, when 
no torque is applied. In figure [8^, we plot the PDF of 59 measured at M ^ together 
with the Gaussian fit of the PDF at equilibrium (continuous line). The two curves 
matches perfectly within experimental accuracy. Thus we conclude that the external 
driving does not perturb the equilibrium distribution of angular fluctuations, so we use: 



The third experimental observation concerns time correlations. In figure [Hb, we plot 
the power spectral density function of 59 when applying an external forcing (o). We 
compare it to the prediction of fluctuation dissipation theorem at equilibrium (Eq. [2]) 
computed using the oscillator parameters. The two spectra are identical, so we can 
confidently use for our system a description in terms of a second order Langevin dynamic 
where the noise term is not perturbed by the presence of the driving. From the power 
spectral density function of 9 (Eq. [2]), we derive the autocorrelation function Rse(r) of 
59 during a time interval r. It is the same at equilibrium and out of equilibrium, and 




(14) 
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Figure 8. a) PDF of the fluctuations 56 = 9 — 6(t) when the torque is applied (o), 
compared with a Gaussian fit of the PDF at equilibrium (continuous line), b) The 
measured spectrum of 69 (o) is compared with the prediction of fluctuation dissipation 
theorem in equilibrium (continuous line). 



decreases exponentially: 

Rse(r) = (58(t + r)59(t)) = exp ( -^P) sin (^| T | + ^ (15) 

where ifj 2 = (ujo) 2 — (l/r a ) 2 and <p is denned by cos(y?) = l/(u r a ) and sin(</?) = ip/u^. 

Thus we observe experimentally that when we drive the system out of equilibrum, 
the angular fluctuations 58 are identical (with respect to the expressions above) to those 
at equilibrium. We verify the same properties for the sinusoidal time prescription of the 
torque, and use equilibrium expression for the correlation function in the next sections. 



7.2. Work distribution 

In Figs. [31 H] and El we see that the PDFs of the work are Gaussian for any integration 
time r and whatever the forcing is. So these distributions are fully characterized by 
their mean value (W T ) and their variance o j1 Vt = = ((W T — (W T )) 2 ). The external 

torque M is determistic, so the mean value of the work done on the system can be 
written as : 



1 



W> = T~r \ M (t')(0(t'))dt'. (16) 
k B T J t . 

We have defined M(t') = M(t') —aM(ti). The value a depends on the time-prescription 
of the torque we apply to the oscillator. Choosing a = 1 describes the linear ramp and 
a = corresponds to the sinusoidal forcing. 
The variance of the PDFs is : 

a 2 Wr = ^^£ l+T £ l+T M(t 1 ).M(t 2 )(^(t 2 )^(t 1 ))dt 1 dt 2 . (17) 
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This expression involves the autocorrelation function of the angular speed 89, 
((S9(t2)59(ti))). Using the expression of the autocorrelation function of angular 
fluctuations (Rse, Eq. (1T51)). we can calculate exactly the expression of (59(t 2 )S9(ti)): 

(59(^)89^)) = - - kBT l . exp (-h^A) sm(V|* 2 -*i| -^)- (18) 



I eB sin(ip) V r « 
We have calculated the mean value and the variance of the PDFs in the three situations 
of interest: stationary and transient cases with a forcing linear in time, and stationary 
case with a forcing sinusoidal in time. Details and results can be found in the appendix. 
In all of the cases, we compare the theoretical PDFs and the symmetry functions with 
the experimental results. We have plotted in Fig.[3],|4]and[6]our theoretical PDFs and the 
corresponding symmetry functions with no adjustable parameters. Within experimental 
error bars, our analytical and experimental results are in excellent agreement. S(w T ) 
is linear in w T because the PDFs of the work are Gaussian. We now want to calculate 
analytically the corrections to the slope T, w (t) for finite time r. For a Gaussian 
distribution, the symmetry function is : 

S(w T ) = ^f— w t = £ic(t) w t . (19) 
a w 

The expression of the slope T, w (t) uses only the mean value and the variance of the 
Gaussian distribution. We define S(r) = (1 — e(r)) -1 , where the correction e(r) is a 
decreasing function of r. We obtain e(r) = for the transient case, which is in agreement 
with a TFT. For the two steady states, there are corrections to the value 1; we find: 



i) Linear forcing 

<r) = — 

ii) Sinusoidal forcing 



(20) 



e(r) = -^ + -^e-^. (21) 

Exact values of the coefficients A, B, D, E, F are given in the appendix. These two 
expressions are in perfect agreement with experimental results as can be seen in Fig. |4] 
and Fig. [6j These corrections depend on the kind of forcing and it is difficult to predict 
their form for an arbitrary forcing. Nevertheless, the two situations we consider are 
useful as building blocks of such an arbitrary forcing, and they provide a very nice test 
of our method. 



8. Heat fluctuations : theoretical predictions 

We now determine an analytical expression of the PDF of the dissipated heat. To do 
so, we make the same hypothesis as in the case of the work (see section 17.11 above), 
and we complete them by additional assumptions to simplify our derivations. We are 
interested by PDFs of the heat for integration time r large compared to r a , so that 
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exponential corrections which are scaling like e~ T l Ta can be neglected. In the case of 
sinusoidal forcing, this is correct after 3 or 4 periods of forcing (r/r Q = 1.64n). Within 
this assumption, 9(U + r) and 9(U) are independent, and so are ^(U + r) and ^-(U)- 
Additionally, as the equation of motion of the oscillator is second order in time, 9 and 
4f are independent at any given time t. We use the technique proposed in [7J . To obtain 
the PDF p(Q T ) of the heat, we define its Fourier transform, the characteristic function, 
as 



P T (s) = / dq T e ls ^p(q T ) (22) 

J — oo 

We then write p(q T ) using Eq. as: 

p{q T ) = J J d9d9P (AU T - Q T , 9(U + r), 9(U), 9(U + r), 0(t 4 )) . (23) 

where P is the joint distribution of the work W T , 9 and ^ at the beginning and at the 
end of the time interval r. This distribution is expected to be Gaussian because W T is 
linear in 9 and additionally 9, 9 and W T are Gaussian. The details of the calculation 
are given in the appendix. 

8. 1 . Linear forcing 

The Fourier transform of the PDF of dissipated heat can be exactly calculated : 

-16 cos(y2) 2 + 4 + 4is(4 cos(y?) 2 + 1] 

+ TT7 2 " 

As far as we know, there is no analytic expression for inverse Fourier transform of this 
function, so for the PDF of dissipated heat. However we can do some comments. This 
expression is very similar to the one found in the case of a Brownian particle [7]. The 
factor (1 + s 2 ) -1 is the Fourier transform of an exponential PDF and this is directly 
connected to the exponential tails of the PDF. Moreover the PDF is not symmetric 
around its mean, because there is a non vanishing third moment. In this expression, 
only two terms depend on r. For large r, this expression reduces to: 

P T (s) = -^exp \-2id 2 — s(l +is)\ . (24) 
1 + s 2 I r a J 

This expression will turn out to be similar to the one obtain with a sinusoidal forcing, 
as we will comment in the next section. 

Both expressions depend on the non-dimensional factor d defined as: 



S^-d 2 is ^2— +is 


2— + 2 







1 M 

d =*/«^t- (25 » 



All moments of the distribution of Q T are linear with d 2 and (9) / y (59 2 ) = d. So d 2 
compares the mean value of the angular speed to the root mean square of the angular 
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speed fluctuations. This coefficient d 2 increases when the system is driven further from 
equilibrium. We consider it as a measure of the distance to equilibrium. In our system 
d is positive, but smaller than 1, so we are out-of-equilibrium but not very far from it 
{d = 0.059). 



8.2. Sinusoidal forcing 

Just like in the experiments, we choose the integration time r to be a multiple of 
the period of the forcing, so (AU T ) = and therefore (W T ) = —(Q T )- Within this 
framework, we find that the PDF of AU T is exponential: 

P(A[/ T ) = ^exp(-|AC/ T |). (26) 

It is independant of r because AU T depends only on 9 and ^ at times U and U + r 
which are uncorrelated. This expression is in perfect agreement with the experimental 
PDFs for all times (see Fig. [7b). Some algebra then yields for the characteristic function 
of Q: 

i(Qr)s - 

The characteristic function of heat fluctuations is therefore the product of the 
characteristic function of an exponential distribution (jr^s) with the one of a Gaussian 

distribution (exp \ i(Q T )s — -f 1 ^ 2 ^). Thus the PDF of heat fluctuations is nothing but 
the convolution of a Gaussian and an exponential PDF, just as if W T and AU T were 
independent. The inverse Fourier transform can be computed exactly: 

1 fal 



P T (s) = — — exp i(Q T )s - -^s 2 (27) 



P(Q T ) = - exp 



e Qr-<Qr) Erfc 




<Qr-(Qr)) ETic *r t \^rrr u w 

where Erfc(x) = 1 — Erf(x) stands for the complementary Erf function. In Fig. [Tfc, we 
have plotted the analytical PDF from Eq. (1251) together with the experimental ones, 
using values of aj^ and (Q T ) from the experiment and no adjustable parameters. The 
agreement is perfect for all values of n, i.e. for any time r. From Eq. ( 1281) . we isolate 
three different regions for S(q T ): 

(I) if Q T > a 2 w + \(Q T )\ = 3\{Q T )\ + 0(1), then S{q T ) = 2 + 0(±). This domain 
of S T corresponds to fluctuations larger then three times the average value. The PDF 
has exponential tails, corresponding to an exponential distribution with a non- vanishing 
mean. 

(II) if Q T < a 2 w - \(Q T )\ = \(Q T )\ + 0(1), then S(q T ) = S(n)g T + 0(±) with 
S 9 (n) = 2 ^9 2 T ^ = E w (n). In this domain, values of the heat are small and heat 
fluctuations behave like work fluctuations. The slope S(r) is the same as the one 
found for work fluctuations. The exact correction to the asymptotic value 1 is plotted 
in Fig. [7^ and again it describes perfectly the experimental behavior. 
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(Ill) for cr^y— \ (Q T ) \ < Q T < <Tyy + \{Q T )\, there is an intermediate region connecting 

2 - 

S(r). 



domains (I) and (II) by a second order polynomial: S(q T ) = 2— ^p-(q T — (1 + wpj)) 2 + 



These three domains offer a perfect description of the three regions observed 
experimentally (Fig. [7H). 

Now, we examine the limit of infinite r in which SSFT is supposed to hold. To 
do so, we distinguish two variables : the heat Q T or the normalized heat q T . Their 
asymptotic behavior are different because the average heat (Q T ) depends on r, more 
precisely it is linear in r. 

We discuss first Q T . The asymptotic shape of the PDF of Q T (Eq. (1281) ) for large 
r is a Gaussian whose variance is er^, the variance of the PDF of W T . Thus, the PDF 
of Q T coincides with the PDF of W T for r strictly infinite. As we have already shown, 
work fluctuations satisfy the conventional SSFT ; therefore heat fluctuations also satisfy 
the conventional SSFT (Eq. (jHJ)). We have found three different regions separated by 
two limit values: the mean and three times the mean. But in the limit of large times 
r, the PDF shrinks and only region (II) is relevant. Region (II) corresponds to small 
fluctuations and it is bounded from above by |(Q r )| with the average (Q T ) being 

linear in r. So all the behavior of the fluctuations of Q T for large r lays in region (II) 
where the symetry function is linear and SSFT holds. 

We turn now to the normalized heat q T . As the average value of Q T is linear in 
r, rescaling by (Q T ) is equivalent to a division by r ; the mean of q T is then 1. This 
normalization makes the two limit values constant. The boundary between regions (II) 
and (III) is 1 + 0{l/r) and the boundary between (III) and (I) is 3 + 0(l/r). The 
function S(q T ) is not linear for large values of q T > 1 but it is linear only in region (II), 
for q T < 1, i.e. for small fluctuations. So SSFT is satisfied only for small fluctuations 
but not for all values of q T , and we obtain for q T a fluctuation relation which prescribes 
a symetry function that is non-linear in q T . 

These two different pictures, in terms of Q T or q T , results from taking two non- 
commutative limits differently. The first description using Q T implies that the limit r 
infinite is taken before the limit of large Q T . The second description does the opposite. 
However, the probability to have large fluctuations decreases with r and experimentally, 
for large r, only the region (II) can be seen, and it is the region in which where SSFT 
holds. 

As we have done in the case of the linear forcing, we introduce a non-dimensional 
factor d such as : 

d = x f^Fl M f" (29) 



\ 




+4(^cob (s ,)] (30) 



The moments of the distribution of Q T are linear with d 2 and, like the linear torque, d 



is equal to the amplitude of 9 divided by y (59 2 ). We consider it also as a measure of 
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the distance to equilibrium. In our system d is positive, but smaller than 1, so we are 
out-of-equilibrium but not very far it : here d — 0.18. 

9. Discussion and conclusion 

In this paper, we have studied the fluctuations of energy input and energy dissipation 
in a harmonic oscillator driven out of equilibrium. This oscillator is very well described 
by a second order Langevin equation. We have performed experiments using a torsion 
pendulum driven out of equilibrium following a stationary protocol in which either the 
torque increases linearly in time, or oscillates at a given frequency. We have also studied 
transient evolutions from the equilibrium state. We have defined the work given to the 
system as the work of the torque applied during a time r. Accordingly we have defined 
the heat dissipated by the pendulum during this time r, by writing the first principle 
of thermodynamics between the two states separated by time r. 

Fluctuations Relations are obtained experimentally for both the work and the heat, 
for the stationary and transient evolutions. 

We have experimentally observed that angle fluctuations of the brownian pendulum 
have the same statistical and dynamical properties at equilibrium and for any non- 
equilibrium driving. From this observation, we have derived expressions for the 
probability density functions of the work and the heat. In our system, fluctuations of the 
angle are Gaussian, and so are fluctuations of the work w r . So the symmetry functions 
S(w T ) of the work are linear, and we have calculated exactly the time-correction to the 
proportionality coefficient between S(w T ) and w T . These corrections match perfectly the 
experimental results, both in the case of a forcing linear in time, and sinusoidal in time. 
We have also computed the analytic expression of the Fourier transform of the PDFs of 
the dissipated heat. For the sinusoidal forcing, we have obtained for the first time an 
analytic expression of the PDF of the heat. This expression is in excellent agreement 
with the experimental measurements. For a torque linear in time, the PDF of the heat 
has no simple expression but its Fourier transform gives insight on the behavior of the 
symmetry function of the heat. It is very similar to the one obtained in the case of a 
first order Langevin dynamics [7J. We emphasize here that our analytical derivations 
are strongly connected to experimental observations on the properties of the noise ; and 
are therefore different from any previous theoretical approach. 

We have introduced a dimensionless variable d which we think is a measure of 
the distance from equilibrium: the average dissipation rate is proportional to d, and 
it increases when the system is further from equilibrium, d is also proportional to the 
strength of the driving and in the Fluctuation Relations, it gives a proper unit to measure 
the amplitude of fluctuations. So d plays the same role as the dissipation coefficient (the 
viscosity in our case) in the Fluctuation Dissipation Theorem at equilibrium. We have an 
expression of d for the two different time-prescriptions we have used. These expressions 
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can be generalized: 



d 2 



(31) 



The numerator corresponds to the solution of the Langevin equation when removing the 
thermal noise term (77 = 0). The denominator corresponds to the variance of thermal 
fluctuations of angular speed 89. 

We thank G. Gallavotti for useful discussions. This work has been partially 
supported by ANR-05-BLAN-0 105-01. 

Work fluctuations 

In this section, we will calculate the mean and the variance of the work given to the 
system in the cases : 

(i) Transient state, linear forcing 

(ii) Steady state, linear forcing 

(iii) Steady state, sinusoidal forcing 



TFT, forcing linear in time 

The torque is : M{t) = M t/r r . The mean value of the angular displacement is the 
solution of Eq. [1] : 



(0) 



M, 



— (e~ t/Ta sm(ipt + 2ip) +ijjt- sin(2</?)) . 



(•1) 



For the work done on the system, the PDFs are Gaussian for all integration time r. The 
mean of the PDF of W T for a given r is: 



(W T ) 



and its variance is: 



Ml 



k B T^Cr? 



-ipr 2 + re r/Ta sm(ifjT + 2ip) 

+ — (e~ T,Ta sin(^r + 3ip) - sin(3y?)) 
^0 



a 



2M G 2 



ksTipCr, 



-^t 2 + re t/t " sin(^r + 2<p) 
+— (e~ T/Ta sin(^r + 3<p) - sin(3^)) 



(•2) 



(.3) 



0*y r = 2{W T ) . 



(•4) 
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SSFT, forcing linear in time 

The torque is : M(t) = M t/r r . The mean value of the angular displacement is the 
solution of Eq. [1] after some r a . Thus the exponential term is vanished: 



(0) 



M. 



— (V>t-sin(2y>)). 



(.5) 



For the work done on the system, the PDFs are Gaussian for all integration time r. The 
mean of the PDF is: 

Ml 



<Wr) 



and the variance is: 



2k B TCr? 



2Ml 



a 



W'r 



-ipr 2 + re r/Ta sin('0r + 2cp) 



H (e T,Ta sin(^r + 3ip) - sin(3</?)) 



From this, we deduce: 



where 



e r 



A = 2 



-L (a - e"^- (a + — 



sin(3<y9) 



V T 

B = 2 sin(V'r + 2y>) , 
£) = 2 sin(^r + 3y>) . 



(.6) 



(•7) 



(.8) 



SSFT, forcing sinusoidal in time 

The torque is M(t) = M Q sin(uJdt). The mean value of the angular displacement is: 

M 



(6) = 6q sm(ujdt + (3) where 6q 



where 



cos(/3) 



1_ I m 



and sin(/3) 



-2 (g) cos(y) 



(.9) 



\ 



[^y^j +4f^cos^ 



;.io) 



For the work done on the system, the PDFs are Gaussian for all integration time r. The 
mean of the PDF is: 

2 

(.11) 



(W n ) = r-=^ -j-r (r/r«) 



k B TC \ p{uj 
and the variance is: 

a 2 n = 2(W n ) + E + Fe~ T/T < 



;.12) 
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where 



(lK)(l + (^) 2 )cos(2/?) 

(co/u o y(r/r a ) { - 1 

[sin(^r + if) cos(2/3) + (uj/uj ) 2 sin(^r - </?) cos(2/5) 



(Lu/uj ) 2 .(r/T a )) 

+ (w/^ )sin(^r)sin(2/3)] . (.14) 



Heat fluctuations 



In this section, we will calculate the Fourier transform of the PDF of the dissipated heat 
in two cases : 

(i) Linear forcing 

(ii) Sinusoidal forcing 

Linear forcing 

We introduce non-dimensional parameters in order to simplify calculations : 



M(t) \ 



* =*/j^e>. (.is) 



The mean value and the variance of x and x can be simply expressed 
d = 



Ml 



CksT uoT r 

(x) = - 2.d. cos(^) , (x) = d , (Sx 2 ) = 1 , (Sx 2 ) = 1. (.16) 
where d is a non-dimensional value. Integrating by part, the work W T can be rewritten: 

W T = d.Lu [(U + t)x(U + r) - Ux(U)] ~ [{U + r) 2 - t 2 ] + W\ 

W* = - (du ) [ l+ x(t')dt'. (.17) 



With these definitions, we obtain : Q T = |A5 T + ^Ax T — W* and (Q T ) = —(W*). Like 
the distribution of W T , the distribution of W* is gaussian for all values of r and we find: 

(W*) = 2d 2 r/r a , 

aly, = 2d 2 (2r/r Q + 1-4 cos(cp) 2 ) . (.18) 

For notational convenience, we introduce a five dimensional vector : Y = (W*,x(ti + 
r),x(ti),x(ti + r),x(ti)). P is Gaussian and is so fully characterized by the covariance 
matrix C defined as: 

= -ft)?) (.19) 
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where denotes the complex conjugate of Z. So the distribution P is written : 

p < y )-\fi5i exp H (y - <F>)Tc " (F - <y>) 



;.2o) 



where Z T denotes the transpose of Z. We suppose that the integration time is larger 
than the relaxation time. Within this assumption, 9{ti + r) and 9{ti) are independent, 
and so are 9(U + t) and 9{tj). As the equation of motion of the oscillator is second order 
in time, 9 and 9 are independent at any given times t. With these hypotheses, we get : 

(5x(U)5x(U + r)) = (5x(U)5x(U)) = (5x(U)5x(U + r)) 

= (6x(ti)Sx(ti + r)) = (6x(ti)6x(ti + r)) 
= (Sx(ti + r)6x(ti + r)) = 0. 

The other coefficients of the covariance matrix are : 

(5W*5x(U)} = (5W*5x(U + r)) = -2dcos(^), 
(5W*5x(U)} = (5W*5x(U + r)> = -d. 

The Fourier transform of the PDF of the heat will be calculated here. We define 
two quantities : 

/ \ 

10 

0-10 

1 

\ -1 



.21) 



.22) 



/ 1 \ 







N 



.23) 



One can write Q T as : Q T = ^Y T NY — e T Y. The Fourier transform can be so written 



PJs) 



M 



dY 



v /(2vr) 5 detC 



exp (M) 



1 



(6Y) T C~ 1 (SY) +is[ -Y 1 NY - e r Y 



We use a new variable defined as : 

Y' = Y - (1 - i sC.iV) _1 ((y) - is.C.e). 

With this definition, the argument in the exponential M can be rewritten : 

1. 



;.24) 



.25) 



M 



7 



-Y' {C 1 - isN)Y' + 7, 



i s 



y [(N(Y) - e) T (l - i sCN)-\(Y) - f sCe) - (Y) T e\ . 



(-26) 



Changing the integration variable to V yields : 

ay 



^ w = /7W exp (-^' r(C_I " isJV)y ') ' exp(7) 



v/(2vr) 5 detC 
exp (7) 
^/det(l -isC.N) 



.27) 
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To make Eq. 1.271 into an explicit expression for P T , the inverse of matrix (1 — isC.N) is 
required in the expression for 7 and its determinant. These are obtained as follows. We 
find: 

/ 1 i s(2dcos((p)) i s(— 2dcos((p)) isd —isd \ 
l-is 

l + is 

l-is 

l+is 

The determinant of the matrix is (1 + s 2 ) 2 . For the inverse of this matrix, we get 

is -d 



isC.N 







/ 



isC.N)- 1 



( 1 




1 

l—i s 














^0 






(2dcosM) ^- s (2dcos(cp)) 



1 

l+is 





1— i s 




1 

l—i s 











d \ 



1+i: 



We now have the material needed to calculate 7. We find : 



.2 r 



7 = -is(W*) 



cx^ + 2c? 2 (l + 4cos(^) 
' As 2 d 2 



Ai sd 2 
1 + s 2 



(4cos((^) 2 - 1) 



+ 



■(4cos(^) 2 + l) 



'.281 



1 + s 2 

The analytic expression of the Fourier transform of the PDF of the heat dissipated 
during a linear forcing is : 

1 



PAs) 



1 



cxp 













-is 















-16 cos(v?) 2 + 4 + Ai s(A cos(y?) 2 + 1] 
1 + s 2 



(-29) 



Sinusoidal forcing 

We will determine in a first time the Gaussian joint distribution P of W T , 9(U), 9(ti + r), 
6(ti) and 9(ti + r). For notational convenience, we introduce a five dimensional vector 
: x = {W T ,6{ti + T),6(U),6(ti + t),6(U)). The PDF P is fully characterized by the 
covariance matrix C. 

C l3 = {{X l -{X l )){X J -{X0) (.30) 

where denotes the complex conjugate of Z. We suppose that the integration time 
is larger than the relaxation time. Within this assumption, 9{ti + r) and 9(ti) are 
independent, and so are 9(ti + r) and 9(ti). As the equation of motion of the oscillator 
is second order in time, 9 and 9 are independent at any given times t. With these 
hypotheses, we get : 

(59(U)59(U + r)> = (6e(U)66{U)) = (59(U)59(U + r)> 

= (59(U)59(U + r)> = (59(U)59(U + r)) 
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= (seiu + ^seiu + r)) = (w T e{u)) 

= (W t 6(U + t)} = (W T 6(U)} 
= (W T 6(U + r)> = 0. 

The covariance matrix is so a diagonal matrix : 



C 



.31) 



/ 
















\ 







k B T/C 





















k B T/C 





















k B T/ I cS 







V 














k B T/ J cff 


/ 



;.32) 



AU T is a function of the positions and velocities at the beginning (tj) and at the end 
(t; + r). Thus, AU T and W T can be considered as independent. The PDF of Q T is so 
the convolution between the PDF of W T which is Gaussian and the PDF of AU T : 

-+oo 

.33) 



/+oo 
Pw T (z)P A U T (Qr + z)dz. 
■oo 



We first calculate exactly the PDF of the variation of internal energy. We have 
shown that 9(ti), 9(ti + r), 9(ti) and 9(ti + r) are independent. The Fourier transform 
of the PDF is : 



where E n 



Pau t (s) = PE p (t l+ T)(s).P Ec{ti+T) (s).P Epiu) (-s).P Ec (u)(- s ) 
1 ;C9 2 and E r 



;.34) 



- 2k^fles9 ■ The distribution of 9 is gaussian with 
variance k B T/C. The distribution of E p and the distribution of E c are the same : 



P E Jx)=P E (x) 



7TX 



exp(— x) 



The Fourier transform of this distribution is : P(s) = (1 



is) 



(.35) 

~ x l 2 . This distribution is 



the same for E p and E c at tj and t{ + r. Thus the Fourier transform of the variation of 



internal energy is 

PauM = 
and the probability is : 

P(AU T ) = 



:i + s 2 



2 6XP( 



IAC/J 



;.36) 

(-37) 



As AU T and W T are independent, the Fourier transform of the dissipated heat can 
be calculated : 



Pr 



exp (i(Q T ) 



2 3 



1 



.38) 



This expression can be inversed because it is simply the convolution between a Gaussian 
distribution and an exponential distribution. So we find Eq. (1281 . 
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